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Abstract: We use the recently obtained two-loop hard thermal loop perturbation the- 
ory thermodynamics functions of a plasma of quarks and gluons to compute the diagonal 
second- and fourth-order quark number susceptibilities. The two-loop hard thermal loop 
perturbation theory thermodynamic functions used are reliable in the limit that the ratio of 
the quark chemical potential to temperature is small. Using this result, we are able to ob- 
tain (semi-)analytic expressions for the quark number susceptibilities at leading- and next- 
to-leading-order in hard thermal loop perturbation theory. We compare the hard thermal 
loop perturbation theory results with perturbative quantum chromodynamics calculations, 
a Polyakov-loop Nambu-Jona-Lasinio calculation, and lattice quantum chromodynamics 
results. Finally, we introduce a prescription for avoiding over-counting in fixed-loop-order 
hard-thermal-loop perturbation theory results for the quark number susceptibilities. 
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1 Introduction 

Dynamical chiral symmetry breaking and confinement are two well-known fundamental 
features of quantum chromodynamics (QCD). Nowadays it is generally believed that at 
high temperature and/or baryonic density strongly interacting matter will undergo chiral 
symmetry restoration and deconfinement phase transitions to a state of matter called the 
quark-gluon plasma (QGP). The study of such phase transitions has become a field of 
both theoretical and experimental interest since information about these phase transitions 
has the potential to provide a more fundamental understanding of QCD itself. In recent 
years, considerable effort has been dedicated to the creation of a QGP in the laboratory. 
The collider experiments currently dedicated to this search are the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Lab (BNL) and the Large Hadron Collider 
(LHC) at the European Organization for Nuclear Research (CERN). Future experiments 
are planned at the Facility for Antiproton and Ion Research (FAIR) at the Gesellschaft fur 
Schwerionenforschung (GSI) facility. In all cases experimentalists collide heavy ions which 
have been accelerated to relativistic speeds in order to create the conditions necessary for 
the creation of a short-lived quark-gluon plasma. 

In the confined/chiral-symmetry-broken phase, quarks are confined inside hadrons 
which possess integer- valued baryon number. In the deconfined/chiral-symmetry-restored 
phase, quarks which have fractional baryon number are free to propagate over larger dis- 
tance scales. This fundamental difference leads to different quark number fluctuations in 
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the two phases [1]. In general, fluctuations of conserved quantities [2], such as baryon num- 
ber, electric charge, strangeness, isospin charge, etc. are considered to be an important 
diagnostic tool for investigating the quark-hadron phase transition in relativistic heavy ion 
collisions [3, 4]. The magnitude of quark number fluctuations can be determined by com- 
puting the quark number susceptibilities (QNS), which measure the response of the quark 
number density to an infinitesimal change in the quark chemical potential in the limit of 
zero chemical potential. In addition, it was also argued recently that the QNS may be 
used to identify the position of the critical end point in the QCD phase diagram [5]. The 
QNS have been extensively studied in the last two decades using a variety of approaches 
including perturbative QCD (pQCD) [3, 6-8], lattice QCD (LQCD) simulations [9-17], 
Nambu-Jona-Lasinio (NJL) models [4, 18, 19], Polyakov Loop-Nambu-Jona-Lasinio (PNJL) 
models [19-22], hard-thermal/dense- loop (HTL/HDL) resummation techniques [23-27], 
rainbow-ladder and beyond-rainbow-ladder approximations of the Dyson- Schwinger equa- 
tions [28] , strong-coupling techniques [29] , functional renormalisation group techniques [30] , 
and quasiparticle models [31]. 

In view of the ongoing experimental and theoretical effort to understand the phase 
structure of QCD, the determination of QNS using a variety of approaches is important. 
Although the study of QCD matter very close to the phase transition requires a nonper- 
turbative description, it is also interesting to explore the behavior of QNS beyond leading 
order (LO) [23-27] within weak coupling expansions which employ state-of-the-art resum- 
mation techniques [32-36]. Unlike LQCD, the weak-coupling expansion can straightfor- 
wardly handle finite density and temperature. As a result, such calculations can provide 
useful information about QNS. Herein we will assume massless quarks. For massless quark 
flavors the QNS are usually defined as 



gr*p 

Xn{T) 



(1.1) 

/u=0 



where V is the pressure of system, fx is the quark chemical potential and T is the tempera- 
ture of the system. Above, n is the order of the QNS and we note that all odd orders vanish 
at fx = due to the charge-parity symmetry of the system. The equation of state (EOS) 
expressed in terms of the pressure V of QCD matter at high temperature and density is an 
important quantity, which we have recently computed to next-to-leading order (NLO) [37] 
at nonzero n and T employing the hard thermal loop perturbation theory (HTLpt) reor- 
ganization of finite temperature QCD [32-36] . In this paper we compute the second- and 
fourth-order QNS from the LO and NLO HTLpt pressure and compare our results with 
conventional perturbative QCD calculations, various LQCD results, and a Polyakov-loop 
Nambu-Jona-Lasinio model. 

The paper is organised as follows: In Sec. 2 we briefly review the HTLpt formalism. 
In Sec. 3 we present (semi-) analytic expressions for the LO and NLO HTLpt pressure. In 
Sec. 4 the second- and fourth-order QNS at LO and NLO are derived. The corresponding 
results are discussed and compared with first principles LQCD calculations, conventional 
pQCD, and a PNJL model. We conclude and give an outlook for future work in Sec. 5. 
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2 Hard Thermal Loop Perturbation Theory 



HTL perturbation theory [34-36] is a reorganization of the perturbation series for hot and 
dense QCD which has the following Lagrangian density 

C = (£qcd + £htl) _ +A£ H tl, (2.1) 

g-^VSg 

where A£htl collects all necessary additional HTL renormalization counterterms and 
£htl i s the HTL effective Lagrangian [32, 33]. The HTL effective Lagrangian can be 
written compactly as 

CuTL = -\(l-6)mlTv(G, a ^jp^^ G A +{ i- 5)im ^^J^ ^ (2 .2) 

where D is the covariant derivative, y = (l,y) is a light like vector and (• • • ) is the average 
over all possible directions, y, of the loop momenta. The HTL effective action is gauge 
invariant, nonlocal, and can generate all HTL n-point functions [32, 33], which satisfy the 
necessary Ward-Takahashi identities by construction. The mass parameters mo and m q 
are the Debye gluon screening and quark masses in a hot and dense medium, respectively, 
which depend on the strong coupling g, the temperature T, and the chemical potential \i. 
In the end we will formally treat the masses mo and m q in (2.2) as being leading-order in 
g in order to make the calculation tractable [34-36]. The n th loop order in the HTLpt loop 
expansion is obtained by expanding the partition function through order 5 n_1 and then 
taking 5 — > 1 [23, 24, 26, 27, 34-38]. In higher order calculations, one usually fixes the 
parameters mo and m q by employing a variational prescription which requires that the 
first derivative of the pressure with respect to both m p and m q vanishes such that the free 
energy is minimized [34, 35, 37]. In the following section, we briefly describe the recently 
computed finite temperature and chemical potential NLO HTLpt pressure [37]. 

3 Leading- and Next-to-leading-order HTLpt Pressure 

Using the above reorganization of finite temperature/density QCD, the two-loop pressure 
for a plasma of quarks and gluons can be obtained in HTLpt by expanding in the ratio 
of chemical potential to temperature through fourth order in fi/T. Below we quote both 
leading-order (LO) and next-to-leading-order (NLO) expressions for the pressure, along 
with the relevant mass gap equations, which are then used to compute NLO QNS at 
various orders. The calculation of the LO and NLO HTLpt pressure at finite temperature 
and chemical potential was presented in an earlier paper [37] and we refer the reader to 
this paper for the detailed calculation. 

3.1 LO Pressure 

The LO HTLpt pressure through 0(g A ) at any \i is [27, 37] 

, tt 2 T 4 f 7d F ( 120 , 240 A 15 - d F , 9 , 9 

V LO = d A — (l + -J- a (l + + _A*j - -ml - 30^ (1 + 12^) ^ 

+ 30^ + 1^1-1+7 + ^^-60^(^-6)^4 . (3.1) 
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Note that apart from the free contribution, no explicit terms proportional to p appear 
through 0(rhq). The dimensionless variables rho, ~>TLq, A, and p are defined as 



m D 



A 



A 



/' 



2vrT ' " q 2-kT 1 " 2vrT ' ^ 2vrT 
At leading order, the weak coupling expressions for the mass parameters are 



2rn2 



m 2 D 



9 Z T 



[c A + s F (l + l2p 2 )] ; 



9 2 T 2 c F 



l+4£ 2 ) 



(3.2) 



(3.3) 



3 l ~" " M ' r,J ' " 9 4 2 
In the expressions for the pressure and masses above we use the standard notation for the 
various Casimir invariants necessary: dp = N c Nf, dA = N 2 — l, sp = Nf/2 , ca = N c and 
cp = (N 2 — 1)/2N C . At LO the variational method for fixing the mass parameters mp> and 



m q can not be used in practice, since it results in only the trivial solution mo 



m n 



0. 



Instead at LO one canonically uses the LO masses listed in Eqs. (3.3) as the lowest order 
"variational solutions" for the mass parameters. 

3.2 NLO HTLpt Pressure and Variational Mass Gap Equations 

The NLO HTLpt pressure through 0[(/i/T) 4 ] is [37] 
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which is accurate up to 0(g 3 ) and nominally accurate to 0(g 5 ) since it was obtained from 
an expansion of two-loop thermodynamic potential in a power series in mjj/T and m q /T 
treating both tud and m q having leading terms proportional to g. Using the result above, 
the mass parameters mo and m q can be obtained using the variational prescription 
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-T'nlo 
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and 
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(3.5) 



which leads to two gap equations for mo and m q , respectively 
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To obtain the second and fourth-order quark number susceptibilities in HTLpt, one 
requires expressions for m£>, -j^mo, m q , and -§^m q at fj, = from Eqs. (3.6) and (3.7). 1 
We list these here for completeness. The result for the limit of the mp gap equation 



Note that odd derivatives with respect to /j, vanish at /i = 0. Fourth-order derivatives at [i = are 
nonzero, however, they appear as multiplicative factors of the gap equations and are therefore not required, 
as we will see below. 
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For m q one obtains 
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In the expressions above, m/j(0) = m/j(T, A, /j = 0), m" D {G) = |^mD(T, A, /x) 
similarly for m g . 
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4 Quark Number Susceptibility 



We are now in a position to obtain the second and fourth-order HTLpt QNS following 
Eq. (1.1). We also note that the pure gluonic loops at any order do not contribute 
to QNS, however, gluons contribute through the dynamical fermions through fermionic 
loops. This makes QNS proportional to only quark degrees of freedom. Below we present 
(semi-) analytic expressions for both LO and NLO QNS. 



4.1 LO HTLpt second-order QNS 

An analytic expression for the LO HTLpt second-order QNS can be obtained using Eq. (3.1) 



o 2 



V lo (T,A,ii) 



dfi 2 

d F T 2 
3 

+ C ^(c A + s F )(log 2 



2 



(2vrT) 2 dp? 



9 



p,=0 
2 



6)(» 



A 



7 ir 2 

- + 7 H 

2 ' 3 



(4.1) 



where the LO Debye and quark masses listed in Eqs. (3.3) and their ix derivatives have 
been used. 



4.2 LO HTLpt fourth-order QNS 

An analytic expression for the LO HTLpt fourth-order QNS can also be obtained using 
Eq. (3.1) 
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(4.2) 



where, once again, the LO Debye and quark masses listed in Eqs. (3.3) and their [i deriva- 
tives have been used. We note that both X2° m (4-1) and x\° i n (4-2) are the same as those 
recently obtained by Andersen et al. [27]; however, the closed-form expressions obtained 
here have not been explicitly listed therein. 
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4.3 NLO HTLpt second-order QNS 

A semi-analytic expression for the NLO HTLpt second-order QNS can be obtained from 
Eq. (3.4) 



X? L °(T) 



8P_ 

d A T 2 



Pnlo(T,A,^) 



1 



(2ttT) 2 dp, 



2 7 7 NLo(r, A,/t) 



A=o 



2 dp a s 

3 d A tt F 



1 + 4 m D (0) + H f In | - 15.662 j (4 - 3C(3)) m 2 D (0) 



+ 44.953 mj(0) + 



t c(3) + (fHy + ln ™ D(0) ' (16 " 12C(3)) 



m£(0) 



(4.3) 



We note that no [i derivatives of the mass parameters appear in (4.3) and, as a result, 
X2 {T) reduces to such a simple and compact form. This is because the second derivatives 
of the mass parameters with respect to fi always appear as multiplicative factors of the gap 
equations (3.8) and (3.9) and hence these contributions vanish. Numerically solving for the 
variational masses using Eq. (3.8) and (3.9) one can directly compute X2 L °(T) from (4.3). 
Alternatively, we have also computed (T) by performing numerical differentiation of 
the pressure in (3.4) which leads to the same result within numerical errors. 



4.4 Fourth-order QNS in NLO 

A semi-analytic expression for the NLO HTLpt fourth-order QNS can also be obtained 
from Eq. (3.4) 
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Figure 1. Left panel: \2 scaled by the free field value for LO (grey band) and NLO (sea green 
band) in 2-loop HTLpt, 4-loop pQCD (sky blue band) [8], LQCD (various symbols) [9, 11, 15], 
and PNJL model (thick purple line) [22] arc plotted as a function of the temperature. The bands 
in HTLpt and pQCD are obtained by varying the MS renormalisation scale (A) around its central 
value by a factor of two. We also used A^g = 268 MeV and 2-loop a s for HTLpt and pQCD. In 
the PNJL model [22] \2 is obtained using an eight-fermion interaction. The Wuppertal-Budapest 
(WB) group [9] data are obtained using the tree-level improved Symanzik action and a stout smeared 
staggered fermionic action with light quark masses ~ 0.035 m s , with m s being the strange quark 
mass near its physical value. The RBC-Biclefeld collaboration [10] used a p4 action whereas the 
MILC collaboration [15] used an asqtad action. In both cases the light quark mass ranges from 
(0.1-0.2) m s . Right panel: Same as left panel but using 3-loop a s and A^jg = 290 MeV. 



where the double derivatives of the mass parameters with respect to fx survive, but the 
fourth derivatives of the mass parameters disappear as discussed earlier. One can now 
directly compute the fourth-order susceptibility by using numerical solutions of the gap 
equations in (3.8) and (3.11). Alternatively, we have also computed Xi LO {T) by performing 
numerical differentiation of the pressure in (3.4) which leads to the same result within 
numerical errors. 



4.5 Results and Discussions 

Computing the different susceptibilities in HTLpt requires a choice of the renormalization 
scale A, the MS momentum scale A^jg, and a specification of the order of the running 
coupling a s used. In what follows we vary the renormalization scale A by a factor of two 
around a central value of A = 2ttT which results in a band that can be used to ascertain 
the level of minimal theoretical uncertainty. The value of A^jg- depends on the order of the 
running coupling chosen and we fix its value from a recent lattice QCD determination [39]. 
The specific value of Aj^g- used in each case is specified in the figure captions. 

In Fig. 1 we have plotted the Nf = 3 second-order QNS scaled by the corresponding 
free gas limit as a function of the temperature. As discussed above, the bands shown for the 
HTLpt and pQCD [8] results indicate the sensitivity of X2 to the choice of the renormali- 
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sation scale A. However, \2 m both HTLpt and pQCD depends only weakly on the chosen 
order of the running of the strong coupling and in turn only depends weakly on A^jg, as 
can be seen clearly from both panels of Fig. 1. The LO HTLpt prediction for X2 seems to 
agree reasonably well with the available Wuppertal-Budapest LQCD data; 2 however, there 
is a sizable variation among different lattice computations [9, 11, 15] considering improved 
lattice actions and a range of quark masses (see caption). However, lowering the quark 
mass (~ 0.035m s , m s is the strange quark mass) nearer to its physical value [9] seems to 
have a very small effect in the temperature range, as seen from the LQCD data. Note that 
for the Wuppertal-Budapest (WB) lattice data shown in Figs. 1 and 2, Ref. [9] provided a 
parameterization of their X2 data 

X2 (T) = e-Cfe/t+W* 8 ) h [tanh(/ 4 1 + / 5 ) + 1] , (4.5) 

where t = T/(200 MeV), h 3 = -0.5022, h 4 = 0.5950, / 3 = 0.1359, f A = 6.3290, and 
/s = —4.8303. The authors of Ref. [9] performed the fit for data in the temperature range 
125 MeV < T < 400 MeV. Using the parameterization above, we have extended their data 
up to 2000 MeV with a step size of 100 MeV in order to allow comparison with the HTLpt 
results at high temperature. The results for \2 obtained using a nonperturbative PNJL 
model [22] which includes an eight-quark interaction are only available very close to the 
phase transition temperature. 

As discussed earlier, purely gluonic loops do not directly contribute to QNS, and 
instead contribute indirectly through dynamical quarks via internal quark loops which 
contribute to the gluonic Debye mass. The LO HTLpt result requires only a resummed 
HTL quark propagator, which in practice corresponds to static external quarks (valence 
quarks) 3 and thus there is no direct dynamical quark contribution in the LO QNS in (4.1) 
except what comes through the LO mass parameter m£> and its derivative in (4.1). We 
note here that x\° m (4-1) overcounts the perturbative order g 2 due to the inclusion of 
the \i derivatives of the LO mass parameters m q and mo which violates the LO mass 
prescriptions as the lowest order variational solutions within one- loop HTLpt. The reason 
is the following: the mass term mf/T 2 (i = q, D) in the one-loop pressure (3.1) contributes 
to LO QNS through the double derivative with respect to [i as ~ [(jn'i) + fmm'l] = 
rriim'l as — > as n —¥ 0. Now, m^m" is nonzero and finite at \l — > 0, which is the 
source of over-counting. As an example, in LO nip the /i dependence appears from a 
dynamical quark loop (internal quark loop in gluon propagator). The second derivative 
of such dynamical quarks at one-loop order leads to (i) a two-loop quark diagram with 
an internal gluon line, and (ii) a mixed two-loop diagram consisting of a pure quark and 
a pure gluon loop. Also, the double derivative of m q with respect to n corresponds to a 
two-loop quark diagram with two internal gluon lines. This would basically lead to higher 
order corrections to LO mass parameters txid and m q which would then spoil the LO 
order mass approximation as the lowest order variational solution. Therefore, taking into 
account the contribution of such two-loop diagrams in LO QNS can lead to an inconsistency 

2 Our result in this case is exactly the same as that obtained recently by Andersen et al [27]. 
3 This is equivalent to the quenched approximation. 
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Figure 2. Same as Fig. 1 except the LO HTLpt QNS is replaced by the prescription specified in 
Eq. (4.6), which avoids the problem of over-counting in the LO HTLpt QNS (see text). 



in the LO mass prescription in one- loop HTLpt and thus overcounts LO QNS. In what 
follows, as an alternative prescription, in the thermodynamic quantities which are related 
through exact thermodynamic relations either by derivatives or integrations with respect 
to fj, or T of a given thermodynamic quantity, one identifies the mass parameters only after 
differentiations or integrations to preserve a given mass prescriptions in HTLpt. 

Following the above prescription we can obtain LO QNS by performing derivatives on 
explicit \x dependence of one-loop pressure in (3.1) as 



(T) 



d F T 2 



1 



3cf 



/ 



2 Vvr 



(4.6) 



where x{ = d F T 2 /3 = N f T 2 . The above expression avoids over-counting in the LO HTLpt 
QNS. We plot the results of this prescription in Fig. 2. Since the over-counting resulted in 
an additional factor of 2 at LO (4.1), Fig. 2 displays only a small change to the LO HTLpt 
QNS compared to Fig. 1, particularly at smaller T. 

The over-counting encountered at LO can be eliminated by pushing the calculation 
to NLO since the derivatives of the NLO mass parameters with respect to fj, vanish (cf, 
Eq.(4.3)) due to the well-defined variational mass gap equations in (3.6) and (3.7). More- 
over, the NLO result also reproduces the correct 0(g 3 ) coefficient which can be seen from 
(4.3) using LO order masses as 



i-i(£ 

2 Vvr 



G 



7T 



(4.7) 



Nevertheless, we see in Fig. 1 and also in Fig. 2 that NLO HTLpt in (4.3) exhibits a 
modest improvement over the pQCD calculation shown, which is accurate to 0(a^lna s ). 
However, the NLO xi is higher than the LO one at higher temperature and it goes beyond 
the free gas value at lower temperatures. It should be mentioned that although the 2-loop 
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T(MeV) T(MeV) 

Figure 3. Left panel: \i scaled by the free field value for LO and NLO HTLpt as given, respectively, 
in (4.2) and (4.4), 4-loop pQCD [8], and LQCD [11, 16] are plotted as a function of the temperature. 
The bands in HTLpt and pQCD are obtained by varying the MS renormalisation scale (A) around 
its central value by a factor of two. We used Aj^g- = 268 McV and 2-loop a s for HTLpt and pQCD. 
Lattice QCD results [11, 16] are represented by symbols. Right panel: Same as left panel but using 
3-loop a s and A^ = 290 McV. 

calculation rectifies the LO over-counting, it does so by pushing the problem to higher order 
in g. The reason can be understood in the following way: in HTLpt the loop and coupling 
expansion are not symmetrical, therefore at a given loop order there are contributions 
from higher orders in coupling. Since the NLO HTL pressure and thus QNS is only strictly 
accurate to order 0(g 3 ) there is over-counting occurring at higher orders in g, namely at 
0(g 4 ) and 0(g 5 ). A next-to-next-to-leading order (NNLO) HTLpt calculation would fix the 
problem through 0(g 5 ) thereby guaranteeing that, when expanded in a strict power series 
in g, the HTLpt result would reproduce the perturbative result order-by-order through 
0(g 5 ). 

In Fig. 3 we plot the fourth-order QNS (xa) scaled by the corresponding free gas value 
for HTLpt as given in (4.2) and (4.4), pQCD, and LQCD. Both the HTLpt and pQCD 
results exhibit a very weak dependence on the choice of order of the running of a s and 
thus A^jg. Nevertheless, the HTLpt results are found to be far below the pQCD result [8] 
which is accurate to 0(agln(a s )) and the LQCD results [11, 16]. Obviously, the LO xa 
in (4.2) suffers from severe over-counting due to the inclusion of \i derivatives of the LO 
mass parameters that spoil the mass prescription in LO HTLpt. Also, the correction to 
Xa when going from LO to NLO is quite large. This is again due to the survival of the 
double derivatives of the mass parameters with respect to fx in (4.4), which again violate 
the NLO mass prescription obtained through the variational mass gap equations in (3.6) 
and (3.7). As discussed earlier if one suppresses these derivatives of mass parameters, one 
would rectify the over-counting in the contributions to xa i n LO HTLpt. With this recipe, 
i.e., performing the explicit \i derivatives of the one-loop pressure in (3.1), the LO xa can 
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Figure 4. Left panel: Xi scaled by the free field value for NLO HTLpt in Eq.(4.9), 4-loop pQCD [8], 
and LQCD [11, 16] are plotted as a function of the temperature. The bands in HTLpt and pQCD 
arc obtained by varying the MS renormalisation scale (A) around its central value by a factor of 
two. We used = 268 McV and 2-loop a s for HTLpt and pQCD. Lattice QCD results [11, 16] 
are represented by symbols. Right panel: Same as left panel but using 3-loop a s and Aj^g = 290 
MeV. 



be obtained 

x\°(T) = xi[l-0(g 6 )], (4.8) 

where the HTL correction appears only in 0(<? 6 ) since there is no explicit /Li 4 dependence in 
the terms up to 0(g 4 ) in one-loop HTL pressure in (3.1) (such contributions only appear 
starting at 0{g & )). 

On the other hand the NLO x& can De obtained through the explicit fx derivatives of 
the NLO pressure in (3.4) as 



Xi 



NLO/ T \ d-A 
{ ' 4vr 2 



d F a s 
s— H s F < 

d A 7T 



12 + 80m?, (0) I In - - 1.0811 I (C(5) - C(3)) 



A 



-6 19.836 +288 In 



31C(5) - 120 




+ lnm D (0) (C(5)-C(3)) 



(4.9) 



where, as expected, the correct order g 2 term appears along with higher orders, viz., 0(g 4 ) 
and 0(g 5 ); however, no 0(g 3 ) term appears. 

In Fig. 4 we show the result of this modified prescription for the NLO X4 (4-9). As 
can be seen from this figure when one avoids the over-counting there is a remarkable im- 
provement compared to Fig. 3 and the result exhibits a better agreement with LQCD 
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data [11, 16]. Since it is only accurate up to 0(g 2 ), higher orders still suffer from over- 
counting. This suggests that one may need to go beyond NLO to NNLO HTLpt in order 
to fix 0(g 3 ), full convergence, and perhaps even better agreement with pQCD. A simi- 
lar situation was seen with the pressure at zero chemical potential for which the NNLO 
HTLpt result provided a dramatic improvement over the NLO result [36]; however, there 
is no guarantee that a similar dramatic improvement will occur for QNS. The only way to 
know is to extend the calculation for the finite temperature and chemical potential HTLpt 
thermodynamic functions to NNLO. Work along these lines is currently underway. 

5 Conclusion and Outlook 

In this paper we have obtained the second- and fourth-order QNS from the NLO HTLpt 
pressure obtained in a high temperature expansion through C[(^i/r) 4 ]. Analytic expres- 
sions were found for both x\ an d Xa within LO HTLpt. Our result for \\° (4-1) is 
in agreement with the results obtained previously by Andersen et al. [27]; however, we 
have found that the perturbative mass prescription used therein slightly over counts the 
LO QNS, since /1 derivatives of the mass parameters appear which should not if a proper 
variational mass prescription is employed. For the same reason the X4 in (4.2) suffers 
from over-counting through all orders to 0(g 4 ). In this work we proposed a prescrip- 
tion which rectifies this over-counting at LO by requiring the necessary derivatives of the 
mass to vanish. The LO result for xi shows reasonable agreement with available LQCD 
data [9, 11, 15]; however, at this point in time there is still a fairly sizable variation of 
this quantity between the different lattice groups. Moving forward it would seem that a 
detailed analysis of the uncertainties in the various LQCD calculations is necessary before 
detailed conclusions can be drawn. 

At NLO we obtained semi-analytic expressions for X2 L ° anc ^ X4 LO and, after numer- 
ically solving the necessary variational gap equations for the mass parameters mo and 
m q , we obtained our results for X^ LO in (4.3) and x^ LO m (4-4). Unlike the LO results, 
our NLO calculation takes into account dynamical quark contributions by including two- 
loop graphs which involve fermion loops; however, they suffer from the same problem that 
NLO HTLpt calculations at zero chemical potential faced: the NLO X2 in (4-3) gets 0(g 3 ) 
correct but the 0(g A ) and 0(g 5 ) contributions are incorrect if they are expanded out in 
a strict power series in g. As a result, our NLO result for X2 scaled to the free limit is 
closer to unity than the corresponding LO result and only shows a weak dependence on the 
chosen value of the renormalization scale. Our NLO result for xa (cf> ec l-(4.4)) in which [i 
derivatives of the variational mass parameters survive is significantly below the pQCD and 
lattice data. As we discussed, this occurs because of severe over-counting since higher order 
corrections contribute through the derivatives of the mass parameters with respect to \x. 
Once the NLO mass prescription is enforced using the method proposed herein, the NLO 
X4 (4.9) becomes accurate to 0(g 2 ), which was missing in LO HTLpt (4.8). Nevertheless, 
no 0(g 3 ) term shows up at NLO in Xi (4-9), whereas higher orders, viz., 0(<? 4 ) and 0(g 5 ) 
suffer from over-counting like NLO %2- 
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As was the case with the pressure at zero chemical potential, fixing this problem will 
require going to NNLO. In the case of the zero chemical potential pressure, performing such 
a calculation resulted in much improved agreement between HTLpt and LQCD calculations 
above ~ 2T C . At the very least a NNLO calculation will fix the over-counting problems 
through 0(g 5 ). Whether going to NNLO will improve the agreement of the HTLpt X2 
and xa predictions with LQCD results will have to remain an open question for the time 
being. Work on the NNLO calculation has begun, but being a NNLO calculation, care and 
patience must be applied in equal measure. 
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